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Abstract. We define a class of algebraic self-adjoint extensions for the Hodge Lapla- 
cian on manifolds with incomplete edge singularities. We construct the associated 
heat kernels on edge manifolds, reviewing the method of signaling solutions by Moo- 
ers [Moo99] for isolated conical singularities. The arguments in the conical case do 
not carry over ad verbatim to the setup of edges. We then establish the heat kernel 
asymptotics for the algebraic extensions of the Hodge operator on edges, and elabo- 
rate on the exotic phenomena in the heat trace asymptotics which appear in case of 
a non-Friedrichs extension. 



1. Introduction 

Unusual new phenomena in the heat trace asymptotics in the setup of singular spaces 
have attracted a considerable interest since the explicit observations by Falomir, Muschi- 
etti, Pisani and Seeley in [FMPS03] as well as by Kirsten, Loya and Park in [KLP08], 
[KLP06] for certain explicit regular- singular operators on a line segment. A general 
discussion of trace expansions for elliptic cone operators has been provided by Gil, 
Krainer and Mendoza in [GKM 10]. 

For the Hodge-Laplace operator on a manifold with an isolated conical singularity, 
new usual phenomena have already been hinted at by Mooers in [Moo99], [Moo96]. 
However, Mooers did not elaborate in detail on the actual heat trace asymptotics, 
but rather observed certain unexpected non-polyhomogeneity properties of the heat 
kernel. The present work is aimed at closing this gap and a derivation of a full heat 
trace asymptotics for certain self-adjoint extensions of the Hodge-Laplacian in a more 
general setup of incomplete edge singularities. 

A complete characterization of self-adjoint extensions for the Laplacian requires a full 
scale elliptic theory of edge degenerate operators, see [Maz91] and [Sch91]. However, 
in this paper we present a class of algebraic boundary conditions, which define self- 
adjoint realizations of the Hodge-Laplacians on edge manifolds, as already employed 
by the author jointly with Bahuaud and Dryden [BDV11] in context of semi-classical 
non-linear parabolic equations on edge manifolds. 

We then review and extend the heat kernel construction method by Mooers [Moo99] 
to the setup of incomplete edge singularities and algebraic boundary conditions, based 
on the construction of the Friedrichs heat kernel on edges by the author jointly with 
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Mazzeo in [MaVe12]. Hereby, we present (simpler) alternative arguments at various 
steps in the construction. 

In this general geometric setup we recover the unusual new phenomena in the heat 
trace asymptotics, observed in [FMPS03], [KLP08], [KLP06] and [GKM10]. It 
should be noted, however, that there the analysis has been performed irrelated to the 
earlier work by Mooers [Moo99], and in the first three instances relies on a very specific 
exact operator structure and Bessel analysis. 

This paper is organized as follows. We first review the basic geometry of incomplete 
edge spaces in §2. We then classify certain algebraic self-adjoint realizations for the 
Hodge Laplacian in §3 and recall from [MaVe12] the asymptotic properties of the heat 
kernel for the Friedrichs self-adjoint extension in §4. We study the signaling problem 
in §5 and §6. The solution to the signaling problem is the central ingredient in the 
construction of the heat kernel for algebraic self-adjoint boundary conditions, which is 
explained in §7 and is basically a revision of [Moo99]. Finally, in §8 we derive the heat 
trace expansion directly from the heat kernel structure. 
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2. Hodge Laplacian on incomplete edge spaces 

o 

Consider a compact stratified space M with a single top-dimensional stratum M m 
and only one other lower dimensional stratum B b , which is therefore a smooth closed 
manifold. By stratification hypothesis there is an open neighbourhood U C M of B 
and a radial function x : U D M — > M, such that U D M is the total space of a smooth 
bundle over B with fibre ^{F) = (0, 1) x F, an open truncated cone over a compact 
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smooth manifold F*. The restriction of x to each fibre is the radial function of that 
cone. 

Resolution of the startum B in M defines a compact manifold M with boundary 
dM, where dM is the total space of a fibration : dM — > B with the fibre F. The 
resolution process is described in detail for instance in [Maz91]. The neighborhood 
U lifts to a collar neighborhood ^ of the boundary, which is a smooth fibration of 
cylinders [0, 1) x F over B with the radial function x. 

Definition 2.1. A Riemannian manifold (M\dM,g) := (M,g) has an incomplete edge 
singularity at B if over the metric g equals g = go + h, where go attains the form 

g \ <%\dM = dx 2 + x 2 g F + (j)*g B , 

where g B is a Riemannian metric on the closed manifold B, g F is a symmetric 2- 
tensor on the fibration dM restricting to a Riemannian metric on each fibre F, and 
\h\ go = 0(x) as x —7- 0. 

Similarly to other discussions in the singular edge setup, see [Alb07], 
[BDVll],[BaVell] and [MaVe12], we consider a slightly restricted class of edge 
metrics and require <fi : (dM, g F + <p*g B ) — > (B, g B ) to be a Riemannian submersion. If 
p G dM, then the tangent bundle T p dM splits into vertical and horizontal subspaces 
as T^dM © T B dM, where T^dM is the tangent space to the fibre of <fi through p and 
T B dM is the orthogonal complement of this subspace. The requirement for to be a 
Riemannian submersion is the condition that the restriction of the tensor g F to T B dM 
vanishes. Moreover, we require that the Laplacians associated to g at each b G B are 
isospectral, and summarize these additional conditions in the definition below. 

Definition 2.2. Let (M, g) be a Riemannian manifold with an edge metric. This metric 
9 — 9o + h is said to be admissible if 

(i) : (dM,g F + <p*g B ) — > (B,g B ) is a Riemannian submersion; 

(ii) the Laplacians associated to g F at each b G B are isospectral; 

(iii) \h\ go = 0(x 2 ) as x — > 0. 

In order to explain the reason behind the admissibility assumptions, let y = 
(yi, yb), b = dimi? be the local coordinates on B lifted to dM and then extended 
inwards. Let z = (z\, zf), f = dimF restrict to local coordinates on F along each 
fibre. Then (x, y, z) are the local coordinates on M near the boundary. Consider the 
normal operator N(x 2 A p ) yo , defined as the limit of x 2 A p on p-forms with respect to 
the local family of dilatations (x, y, z) — >■ (Xx, X(y — y ), z) as A — > oo. Under the first 
admissibility assumption, N(x 2 A p ) yo is naturally identified with s 2 times the Hodge 
Laplacian on p-forms on the model edge x F x with incomplete edge metric 
g ie = ds 2 + s 2 g F + \du\ 2 . 

The second condition on isospectrality is posed to ensure polyhomogeneity of the 
associated heat kernels when lifted to the corresponding parabolic blowup space. More 
precisely we only need that the eigenvalues of the Laplacians on fibres are constant 
in a fixed range [0,1], though we still make the stronger assumption for a clear and 
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convenient representation. The final condition is of technical origin, placed in several 
other instances, compare [BaVell] and [MaVe12]. 

At any hypersurface S a = {s = a} of the model edge, TS a = T{R b x F) splits into 
the sum of a 'vertical' and 'horizontal' subspace, where the first is the tangent space 
to the F factor and the second is the tangent space to the Euclidean factor M. b . This 
splitting is orthogonal, and induces a bigrading 

(2.1) A p (S a ) = A j (R b ) <g> A\F) := M>\S a ). 

j+i=p j+i=p 

Let n j ' l (S) denote the space of sections of the corresponding summand in this bundle 
decomposition. We write the normal operator N(x 2 A p ) yo with respect to a rescaling 
of the form bundles, employed also in [BrSe87]. More precisely, for each j,l with 
j + I = p, we define 

: C^(R + ,Q j ' l -\S) © n j '\S)) -> %(R b x tf(F)), 

( V , fi) !—)• s 1 - 1 -}' 2 ^ Ads + s l ~ f/2 ^ 

where the lower index refers to the compact support of functions and differential forms, 
away from {x = 0}. We denote by $ p the sum of these maps over all j + 1 = p. Like in 
the special case of conical singularities, the resulting transformation 

$ P :L 2 ([0,1],L 2 ( n j > l -\S) ®Q j > l (S),K(0) + \du\ 2 ),ds)-+ L 2 (n p (R b x tf(F)),0ie), 
j+l=p 

is an isometry, and a calculation yields 

(2-2) % 1 [s" 2 iV(* 2 A p )] $ p = (-^ + I(X - 1/4)) + A R>> 

where A is the nonnegative self-adjoint operator on fi /_1 (F) © f^(F) given by 

A HF + (I-(/ + 3)/2) 2 2{-l) l S l>F 

2(-l) l d l . hF A,, F + (/-(/-l)/2) 2 

Under this transformation the indicial roots of A* have a particularly simple form. 
Writing the eigenvalues of A as z/ 2 , Uj > 0, with corresponding eigenform 0j, the corre- 
sponding indicial roots of A p are given by 

(2-4) 7+ = ^. + !, 7 - = - u . + i. 

A similar rescaling $ p using powers of the defining function x makes sense in each 
local coordinate chart near the boundary DM. Rescalings on different local coordinate 
charts are equivalent up to a diffeomorphism. 

Conjugating by <3> p as before, the rescaled Hodge-Laplacian A* is a perturbation of 
(5.4) with higher order correction terms determined by the curvature of the Riemannian 
submersion : dM — > B and the second fundamental forms of the fibres F. By an 
abuse of notation we denote the rescaled operator by A p again, if there is no danger of 
confusion. 



(2.3) A 
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3. Algebraic boundary conditions on edges 

The exposition of this section has been presented in the recent work [BDV11] of 
the author jointly with Bahuaud and Dryden, and inserted here for convenience of the 
reader and notation fixing. The Hodge Laplacian of an admissible incomplete edge 
space (M,g) is need not be essentially self-adjoint on its core domain C^Q P (M), and 
boundary conditions at the edge are posed to obtain self-adjoint extensions. 

For spaces with isolated conic singularities, this was first studied by Cheeger [Che83] 
using the language of ideal boundary conditions. Further and more systematic studies 
were presented by Bruning-Lesch [BrLe87], see also [Moo99] and [KLP08] with an 
appendix by the author, for a full classification of self-adjoint boundary conditions in 
the conical setting and for results about the associated heat equation. 

In case of isolated conical singularities the extension problem is finite dimensional. 
When the edge B has positive dimension, the requisite analysis is more intricate and 
here we can specify only a class of algebraic boundary conditions for the Hodge Lapla- 
cian. Let A* = A p denote the rescaled Hodge Laplace operator acting on differential 
forms of degree p on the incomplete edge space (M, g) with an admissible incomplete 
edge metric g. Consider the space of square-integrable forms L 2 Q P (M), with respect to 
g. The maximal and minimal closed extensions of A p are defined by the domains 

^ max (A p ) := {u G L 2 n p (M) | A p u G L 2 n p (M)} } 
(3.1) ^min(Ap) := {u G ^ max (A p ) | 3uj G C™VL P such that 

Uj — > u and A p Uj — > A p u both in L 2 Q P }. 

where A p u G L 2 is initially understood in the distributional sense. 

Lemma 3.1 ([MaVe12]). Let (M,g) be an incomplete edge space with an admissible 
edge metric. Consider the increasing sequence of eigenvalues v 2 G [0, 1), j = 1, • • • ,p, 
of the tangential operator A in (2.3), counted with their multiplicities. The associated 
indicial roots of A p are given by 7^ = ±Uj + 1/2. Any u G ^ max (A p ) then admits a 
weak asymptotic expansion as x — > 

g 

i=i 

where the leading order term of each ip^ is the corresponding solution of the indicial 
operator. More precisely, let <pj denote the normalized v 2 -eigenforms of the tangential 
operator A at y. Then 

+ ( VZQogxWjiziy), Vj = 0, 

Wj{x,z;y) = x 1 ! : <j>j{z;y), (x, z;y) = < 

I x 1 ! (1 + ajX)(j)j(z; y), Vj > 0, 

with a,j G M uniquely determined by A p . The coefficients dj[u] are of negative regularity 
in y and the asymptotic expansion holds only in a weak sense, i.e. there is an expansion 
of the pairing f B (u(x, y, z)x{y)) g B dy for any test function x G Q*(B). 
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The set of all closed extensions of A p is in bijective correspondence with the closed 
subspaces of the quotient ^ max / ^mi n ; furthermore, since A p is symmetric on the core 
domain C^°Q P (M), self-adjoint extensions are in bijective correspondence with the La- 
grangian subspaces of this quotient with respect to certain natural symplectic form 
induced from the boundary contributions in an integration by parts formula. However, 
in contrary to the case of isolated conical singularities, the expansion in Lemma 3.1 
holds only in the weak sense, so that self-adjointness of such domains follows only by 
the a mollification argument, as we provide below. 

We proceed with the definition of algebraic boundary conditions which define a sub- 
class of self-adjoint extensions of A p by specifying algebraic relations between the co- 
efficients Cj . Note that though we performed the discussion above under rescalings $ p 
in each local coordinate neighborhood, the partial weak asymptotic expansion of the 
rescaled solutions in ^ max (A) is invariant under coordinate changes x'(x, y, z),y'(x, y, z) 
and z'(x,y,z), with x f (0,y,z) = and y' lifted from the base, so that y'(0,y,z) is in- 
dependent of z. Thus, any specification of algebraic relations between the coefficients 
and hence also the characterization of algebraic boundary conditions for A p will be 
globally well-defined. 

Lemma 3.1 asserts that the quotient A q = ^ max /-^min can be identified with the 
2g-dimensional vector space spanned by solutions {^}j =1 . The symplectic form co q on 
Aq is defined by co q (u,v) := (u,A p v)l2 — (A p u,v)l2 on polyhomogeneous u,v G SUax, 
and explicitly amounts to 

/ 32) u^t^j) = -w q y>j^t) = = 

Following ([Moo96], Section 7), consider now a q x g-matrix T = (Ty) G Matr(g, A g ) 
with diagonal entries given by Tjj = bjjipj + 0jj4>j~, and the off-diagonal entries = 
9ijipj~. The coefficients G R are such that either bu = 1, or bu = 0, where in 

the latter case we require 9a = 1 and 9^ = for i ^ j. We refer to such T = (Ty) G 
Matr(g, A q ) as the Lagrangian matrix. 

Definition 3.2. For any Lagrangian matrix T = (Ty) G Matr(g, A q ) we define the 
associated algebraic domain of the Hodge Laplacian A p by 

^r(A p ) := {u G ^ max (A p ) | V i=lj .. i9 : ^ c^(c+[u]^+ + cJ^^Ty) = 0}. 

i=i 

The algebraic boundary conditions provide a full characterization of self-adjoint ex- 
tensions of the Hodge Laplacian on cones, cf. ([Moo96], Section 7) and also [KLP08]. 
In the setup of incomplete edge, algebraic boundary conditions define only a subclass 
of self-adjoint extensions. A full classification of self-adjoint extensions of the Hodge 
Laplacian on incomplete edges requires a detailed analysis of the elliptic theory of edge 
differential operators and is out of the scope of the present discussion. 



zz/j, Vj > U, 

1, Pa = 0, 
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Below we provide a proof of self-adjointness of @r(A p ). Self-adjoint extensions of A p 
are in one to one correspondence with the Lagrangian subspaces of (A q ,u q ). However, 
standard arguments from the conical setup do not apply here directly, since due to 
weakness of the asymptotic expansion in Lemma 3.1 the symplectic form u q can be 
evaluated explicitly only on polyhomogeneous u, v G £t max and not on the full A q . Hence, 
even symmetry of £^r(A p ) is not obvious here and requires a mollification argument. 

Proposition 3.3. f^r(A p ) defines a self-adjoint extension of A p . 

Proof. We first prove that A p is symmetric on 5?r(A p ). Consider any solution w G 
3>t{A p ). Let be a cut-off function supported in a local coordinate neighborhood 
(x,y,z). Then, still u := w ■ (ft G S>y{A p ). Our first observation is that these locally 
supported elements u G ^r(A p ) can be approximated by f^r(A p ) fl ^hg-elements in the 
graph norm. Consider a coefficient Uj of the form- valued u, a test function if) G C°°{B) 
and write 

(ui*if>)(x,y,z)= / ui(x,y - y, z)ip(y)dy. 

J B 

We can assemble functions uj *ip locally back into a differential form, which we denote 
by u * ip. The convolution u * ip inherits the expansion as x — > from u, and hence 
u * ip G @r(A p ). Moreover, due to pairing with ip G C°°(B), the coefficients cf[u * ip] 
are now smooth in y and hence u* ip G 5?r(A p ) n £^ p hg- Specify now ip to be a bump 
function, compactly supported around the coordinate origin of (y) with ip(0) = 1, and 
consider a sequence ip e (y) := e~ b ip(y/e). Then with respect to || ■ 

^ e (.) = •) ^(0) = 1, as e 0. 

Set u e = u * ip e . Observe for each coefficient u e j = ui * ip e 

\\u t ,i - m 7 || L 2 = \\ui (i>(e-) - || L 2 < HV^(e-) - 1 ||oo ||«/||l2 ->• 0, as e -)• 0. 

This proves u e — > u in L 2 as e — >■ 0. Moreover, A p w e = (A p u) * and hence by 
exactly the same argument as above, A p u e — > A p u in L 2 as e — > 0. Thus, we can indeed 
approximate any locally supported u G 3>t{A p ) by a sequence (u e ) C 2>t(A p ) fl^hg in 
the graph norm. 

By construction of the Lagrangian matrix T, the Hodge Laplacian A p is symmetric 
on ^r(Ap) n ^phg by exactly the same arguments as in the case of an isolated conical 
singularity, see ([Moo96], Theorem 7.6) and also [KLP08]. This is due to the fact 
that polyhomogeneous elements admit an asymptotic expansion in the strong sense 
and integration by parts arguments apply. 

Now symmetry of A p on i^r(A p ) is obtained as follows. Consider a partition of unity 
{<t>a)aeA subordinate to coordinate charts of M. For any /, g G f^r(A p ) we can write 

(A p /, g) L 2 - (/, A p g) L 2 = (( A p/ ; 9 ■ 0«)l 2 - (/, \(g ■ 4>o))i?) = 0, 

where the last equality follows by approximating each / and g ■ <p a locally by elements 
of ^r(Ap) fl ^4>h g in the graph norm, and A p is symmetric on ^ r (A p ) fl ^hg- 
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In order to deduce self-adjointness on 5?r(A p ) it now suffices to show 

^(A; r ) = {fe ^ max (A p ) | Wg G ^ r (A p ) : (A p f,g) L2 = (f,A p g) L i} C ^ r (A p ). 

Let / G £^(A* r ). Then in any local coordinate chart we may consider a locally sup- 
ported g G ^r(Ap) fl =e4>hg, with its coefficients cf[g] being smooth with compact sup- 
port in M. b . Then repeating ([Moo96], Proposition 7.4) ad verbatim, we deduce from 
(A p f,g) L 2 = (f,A p g) L 2 that the coefficients c„- [/] in the weak expansion of / in that 
coordinate neighborhood, satisfy the algebraic conditions of @r(A p ). Consequently, 
/G^ r (A p ). □ 



4. ASYMPTOTICS OF THE HEAT KERNEL ON EDGE MANIFOLDS 

In ([MaVe12], Proposition 2.5) jointly with Rafe Mazzeo we have identified the 
Friedrichs extension of A p as the algebraic self-adjoint extension associated to V — 
diag(/0i", .., ipq). We denote the Friedrichs extension of the Hodge Laplacian by A p 
and explain here the polyhomogeneity properties of its heat kernel H^ near the edge. 

We begin by recalling the definition of conormal and polyhomogeneous distributions 
on a manifold with corners, see [Mel93] and [Mel92]. 

Definition 4.1. Let 2U be a manifold with corners, with embedded boundary faces, and 
{(Hi, Pi)}f =l an enumeration of its boundaries and the corresponding defining functions. 
For any multi-index b = (6i,...,6jv) £ ^ we write p b = p\] . . . p$ . Denote by 
Vb(QH) the space of smooth vector fields on 2U which are tangent to all boundary faces. 
A distribution w on Q2J is said to be conormal if w G p b L°°(W) for some b G 
and V 1 ...V l w G p b L°°(W) for all V, G V b (W) and for every I > 0. An index set 
£i = {(7,p)}cCxN satisfies the following hypotheses: 

(i) Re (7) accumulates only at +00; 

(ii) for each 7 there exists P 7 G No such that p < P 7 for each (j,p) G Ei\ 

(iii) if (7,p) G E i: then (7 + j,p') G Ei for all j G No and < p' < p. 

An index family E = (Ei,...,En) is an A^-tuple of index sets. Finally, we 
say that a conormal distribution w is polyhomogeneous on 2U with index family 
E, denoted w G =e^J; g (QU), if w is conormal and if near each Hi we have w ~ 
J2(-y P )eEi a 7,pP7(l°SP«) P ' as Pi ~~ ^ 0' with coefficients a 7iP conormal on H { and poly- 
homogeneous with index Ej at any Hi fl Hj . 

Consider a local coordinate chart (x, y, z) in the collar neighborhood f C M of 
the boundary dM and note that as in (2.1), the tangent bundle of the hypersurface 
%x = {x = xq} fl tft splits into the sum of a 'vertical' and 'horizontal' subspace, where 
the first is the tangent space to the fiber F 2 -factor and the second is the tangent space 
to the base .By-factor. This splitting is orthogonal, and induces a bigrading 



(4.1) 



= k\B y )®k\E z ) := A l '\W). 

l+k=p l+k=p 
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Then, locally, under the rescaling transformation $ p , the heat kernel is a distribu- 
tion on M\ := K + x M 2 , acting over M. + x ^ 2 on and taking values in the sections 

r((0,l), (A''^ 1 (^)©A / ' fc (^)) 

\ l+k=p 

Consider local coordinates (t, (x,y,z), (x,y,T)), where (x, y,z) and (x,y,z) are coordi- 
nates on the two copies of M near the boundary. The heat kernel has non-uniform 
behaviour at the submanifolds 

P = {(t, (x, y, z), (x, y, z)) e Ml \ t = 0, x = x = 0, y = y}, 

D = {(t, (x, V, z), (x, y, z)) e Ml 1 1 = 0, x = x, y = y, z = z}. 

The asymptotic behaviour of near these submanifolds of M 2 depends on the angle 
of approach to these submanifolds. This dependence on the angle in the asymptotics 
of the heat kernel is conveniently treated by introducing polar coordinates around P 
and D. Geometrically this corresponds to appropriate blowups of the heat space Ml, 
such that the corresponding heat kernel lifts to a polyhomogeneous distribution on that 
blowup space, in the sense of Definition 4.1. 

To obtain the correct blowup of Ml in the case of an incomplete edge metric, one first 
blows up of the submanifold P parabolically. The notion parabolic refers here basically 
to the fact that we treat yt instead of t as a smooth variable. 

The resulting heat space [Ml,P] is defined as the union of M 2 \P with the interior 
spherical normal bundle of P in Ml . The blowup [Ml , P] is endowed with the unique 
minimal differential structure with respect to which smooth functions in the interior of 
M 2 and polar coordinates on Ml around P are smooth. This blowup introduces a new 
boundary hypersurface, which we denote by ff (the front face). 

The actual heat-space blowup ^£ is obtained by an additional parabolic blowup of 
[Ml, P] along the diagonal D lifted to a submanifold of [Ml, P}. The resulting blowup 
is defined as before by cutting out the submanifold and replacing it with its spherical 
normal bundle, which leads to a new boundary hypersurface, the temporal diagonal (td). 
The blowup ^£ = [[Ml, P], D] is a manifold with boundaries and corners as depicted 
in Figure 1. In addition to the new boundary faces ff and td, there are three other 
boundary faces rf, If, tf which arise from the lifts of {x = 0}, {x = 0}, {t = 0} C M 2 , 
respectively. 

In Figure 1 we have also written out some of the appropriate projective coordinates 
on ~-# 2 . Near the top corner of ff away from tf the projective coordinates are given by 

(4.2) P = Vt, £ = -, £ = -, u = , y, z, z, 

P P P 

where in these coordinates p, £, £ are the defining functions of the faces ff, rf and If, 
respectively. For the bottom corner of ff near If, the projective coordinates are given by 

t x y — y _ 

(4.3) r = —, s = -, u = , x, y, z, z, 

X z X X 
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Figure 1 . Heat-space blowup for incomplete edge metrics 



where in these coordinates T ; 5, X £1X6 the defining functions of tf, If and ff, respec- 
tively. For the bottom corner of ff near rf the projective coordinates are obtained by 
interchanging the roles of x and x and are given by 

/ a a\ t x y ~y ~ ~ ~ 

(4.4) T = ^ s = ~> u = x, y, z, z. 

JU %Kj JU 

The projective coordinates on near the top of td away from tf are given by 

(4.5) i] = y/f, S = - — -, U = -, Z = - — -, x, y, z. 

7] 7] 7] 

In these coordinates tf is the face in the limit \(S, U, Z) \ — > oo, and ff and td are defined 
by x and rj, respectively. The blowup heat space is related to the original heat 
space via the obvious 'blow-down map' (3 : ^£ —> which in local coordinates 
is simply the coordinate change back to (t, (x, y, z), (x, y,z)). 

We can now state the asymptotic properties of the (rescaled) Hjt as a polyhomoge- 
neous distribution on the blowup studied by the author jointly with Mazzeo in 
[MaVe12]. 

Theorem 4.2. ( / [MaVe12] ; Theorem 1.2) The heat kernel lifts under the rescal- 
ing $ via the blowdown map j3 to a polyhomogeneous distribution (3*H^ on with 
asymptotic expansion of leading order (—1 — dimi?) at the front face ff and (— dimM) 
at the diagonal face td, with index sets at the right and left boundary faces given by the 
indicial roots 7 > 1/2 (see (2 A)) of the Hodge Laplacian. 

In fact, [MaVe12] went beyond the heat kernel construction, establishing an ana- 
logue of the even-odd calculus for edges with consequences for metric invariance of 
analytic torsion. Below we require a rather detailed understanding of the heat kernel 
asymptotics and are led to provide a short overview of the heat kernel construction in 
[MaVe12] for the Friedrichs extension. 

Definition 4.3. Let S = (E^, E r f) be an index family for the two side faces of We 
define WfX{M) to be the space of all (rescaled, i.e. considered under the rescaling $) 

operators P with Schwartz kernels Kp which lift to polyhomogeneous functions Kp on 
, with index family {(-b - 3 + I + j, 0) : j G N } at ff, {(-m + p + j, 0) : j G N } 
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at td, at tf and £ for the two side faces If and rf of When p = oo, E t d = 0- For 
simplicity we usually denote the lifted Schwartz kernel simply by Kp again. 

As the name suggests, kernels in the calculus V&e-'h (M) may be composed and we 
state the corresponding composition result from [MaVe12]. 

Theorem 4.4. For index sets E^ and E' T{ such that E^ + E' T{ > —1, we have 

where the index sets at the side faces of ^M 2 amount to 

Pi f = E[ { U (E l{ + I') U {{z,p + q + 1) : 3 (z,p) G E[ h and (z, q) G {E xt + /')}, 
P rf = E ri U (E' T{ + I) U {{z,p + q + 1) : 3 {z,p) G £ rf , and (*, g) G + /)}. 

We continue under the rescaling introduced in §2 and do not explify the transfor- 
mation in the notation below. 

The heat kernel is a solution operator to C = d t + A. We lift tC to in projective 
coordinates (4.5) and (4.4) and obtain for its restrictions to the temporal diagonal and 
the front face 

N tA (tC) = ~r)d v + A s>u , z - ^(Sd s + Ud v + Zd z ), 
N e (tC) = r(d T - d 2 + s-\A - 1/4)) + Af ) , 

where A is the rescaled operator (2.3) acting in the variable z. Note that Ng(tC) acts 
tangentially to the fibres S^ 1 of ff, i.e. it involve no derivatives with respect to (x, y,z). 



In order to define a parametrix which solves the heat equation to first order at the 
front face, note that for any H G ^'^(M) with Ng(H) as its leading coefficient in the 
front face expansion, we have 

(4.6) N s (tC oH) = N s (tC) o N S (H). 

To make this vanish, we choose Ng(H) to equal the fundamental solution for the heat 
operator Ns(t£), which is simply the heat operator for the Friedrichs extension of the 
Hodge Laplacian (5.4) on the product (R+ x F x R^, ds 2 + s 2 gp + du 2 ). In other words, 
we set 

N S {H) = H v {t,s,z,s = l,z)H* b (r,u,u = 0), 

SS , / SS \ I S 2 + S 2 



H w (t } s } z } s,I} = }_^ —I v {-) exp ^ —J <f> v (z) ® cf>t(z). 

where we sum over eigenvalues v 2 of the tangential operator A in (2.3) and <p u denote 
the corresponding normalized eigenforms. I u denotes the modified Bessel function of 
first kind. The factor is the Friedrichs heat kernel of the rescaled operator 
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as in [Che83], [Moo99] and [Les97]. The expansion of N S (H) near td is of order 
(—to), where m = dimM. The leading order of the heat kernel H at ff does not follow 
from (4.6) but rather from the initial condition that the heat kernel reduces to the 
delta function as t = 0. Because of the homogeneity of the delta function, the leading 
exponent of the heat kernel at td is (—to). Similar considerations hold at the front face, 
where taking into account the incomplete edge volume form the leading exponent at ff 
is (-6-1). 

Thus we may define an initial heat parametrix G ^l^f(M) by setting it to 
equal p^N^II) near the front face and extending it smoothly to zero the interior of 
The expansion of near td has a universal (euclidean) form and the error of 
the initial parametrix at the temporal diagonal td is solved away following [Mel93], 
uniformly in pg. This defines G (-^0 solving the heat equation to first order 
at the front face and infinite order at the temporal diagonal td. In view of the explicit 
formula for Ng(H) and this proves the following 

Proposition 4.5. There exists an element H^> G *ff^\{M), where £ = (E\f,E r f) and 
the index sets En = E T f are given by {{y + 1/2 + k, 0) | v 2 G spec(A), k G N}. Moreover 
tCH^ = P« G ^1°^ (1) (M) ; where = (%, E Tt - 1), and 

\imH^\t, x, y, z, x,y,z) = 5(x — x)5(y — y)5(z — z). 

i-s>0 

In the next construction step we choose a slightly finer parametrix with an error 
which vanishes to infinite order along rf as well. 

Proposition 4.6. There exists an element J G ^g-if (M), where £' = (Eu, E r { + 1) 
and ■= HW + J is such that t£H® = G ^'^(Af) and lim^ = Id. 
The identity operator Id corresponds to the kernel S(x — x)5(y — y)S(z — z). 

Proof. This step proceeds exactly as in [Moo99]. The error term P^ from the previous 
step has an expansion along rf, and in order to eliminate a term s 7 a (we are using 
projective coordinates (4.4)) its asymptotic expansion, it is only necessary to solve the 
indicial equation 

(-d 2 s + s~ 2 (A - l/4))u = {r~ l a) . 

This is because all other terms in the expansion of tC at rf lower the exponent in s by 
at most one, while the indicial part lowers exponent by two. Note that r, u, x, y and z 
only enter this equation as parameters. We may solve this equation using the Mellin 
transform. The solution is polyhomogeneous in all variables, including the parameters 
and has leading order 7 + 2 at rf. 

Let J denote a kernel which has asymptotic sum at the right face with terms equal to 
the ones obtained iteratively as solutions to to the indicial equation above. Since r and x 
enter the indicial equation as parameters, we can assume that J vanishes to first order at 
ff and to infinite order at tf, as does the error 

p( i) e $ 3,^'( M ) )kt has 

an expansion 

at rf of two orders higher. So indeed J G ^f(M). We see that = + J has 
all the desired properties, with a new error term P^ G \^^^ El{ '°° (^M) . □ 
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Consider the kernels as convolution operators in time. Then, our parametrix solves 
CH^ 2 ' = Id J t-t~ 1 P ( - 2 \ The final stage in the parametrix construction is then to consider 
the formal Neumann series 

oo 

(id + r i pwy i = id + Y^i-t^p^y := H + po) ) 

3=1 

where t~ l P^ e tf^^- 00 and Theorem 4.4 yields {t~ l P^Y e ^^'*' °. In other 
words, successive composition of t~ 1 P ( - 2 ^ with itself produces an operator which vanishes 
to higher and higher order on ff. In fact, a slightly finer analysis, see [Mel93] shows that 
this formal series is actually convergent. We refer to these sources for the necessary 
estimates, both in the case of compact manifolds and for manifolds with cylindrical 
ends, but it is easily seen that everything there transports immediately to this setting 
(and many other ones as well). Indeed, this is a general feature of such Volterra series. 

By the arguments in [MaVe12] the exact heat kernel is then given by 

oo 

= H {2 \ld + P^) = - * r l P® + * J2(- rlp{2) ) j 

3=2 

= + (J + * + 0{p\ +h ), Ps -> 0. 

Recall the admissibility assumption in Definition 2.2 (iii), which requires g = go + h 
with \h\ go = 0(x 2 ) as x —> 0. We can therefore separate £ into the leading order term 
Nff(C), second order term £', comprised of derivatives {d y ,d y d z } which arise from the 
curvature of the fibration : (dM,g F + (p*g B ) —> (B,g B ), and the higher order terms 
C", which arise from h and do not lower the front face asymptotics. Consequently 

(4.8) Hjf = + (J + * + 0(p^ 6 ), p ff -)• 0. 

5. Solution to the model signaling problem 

We now proceed with the first step in the construction of the heat kernel ifp for an 
algebraic self-adjoint extension Ap. The fundamental idea is to add further terms to 
the heat kernel for the Friedrichs extension of the Hodge Laplacian, which correct 
the asymptotic behaviour of the kernel at rf and If to satisfy the boundary conditions 
of f^(Ar). These additional terms are obtained from the signaling solution, which 
is explained and constructed out of ifjr below in §6. The present section provides a 
preliminary discussion of the signaling problem for a prototype of a model edge l v @ A R b, 
where l u := —d 2 + x~ 2 {y 2 — 1/4), v 6 [0,1) is a regular- singular differential operator 
acting on C^°(0, oo). 

The maximal domain 3? max (l v ) and the minimal domain ^ m in(L) for lv are defined 
parallel to (3.1), and provide the maximal and minimal closed extensions of the regular 
singular operator l v in L 2 (0, oo). As a special case of Lemma 3.1, we refer for instance 
to ([KLP08], Proposition 3.1) for an explicit argument, any u G @ max .(l u ) admits a 
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partial asymptotic expansion 

u ~ c + [u] 4>„(x) + c~[u] 4>~j(x) +u, as x — > 0, 

(x~ u+1/2 v G CO 1) 
0+( x ) = x ^l\ 0;( x ) = ' m£ SvUQ. 

[y/x\0g[X), V = 0, 

The weak expansion of solutions in ^ max (/ jy © A R b) is a parallel structure. The signaling 
solution u(-,t) G 3! ms Jl v © A R &),t G [0, oo), is defined for any given h G L X (R + x R 6 ) 
which is -^io C [0, oo) fl C°°(0, oo) in the time variable, as a solution to the so-called 
signaling problem 

(dt + l v ®Agb)u(x,y,t) = 0, u(x,y,0) = 0, 
^ ' ' c~(u(-,t)) = h(t), t > 0. 

Note that u(-,t) cannot take values in a fixed self-adjoint extension of l v © A R t, since 
by uniqueness of solutions to the heat equation, u(x, y, 0) = then implies u = 0. The 
signaling solution is in fact also unique, since {u G St max (l u © A R t) | c~[u] = 0} defines 
the Friedrichs self-adjoint extension of l v © A R t>, which we denote by Lf © A R 6. 

The heat kernel of the Friedrichs self-adjoint extension Lf of the regular singular op- 
erator l v and its restriction to {x = 0} are explicitly given by, see ([Les97], Proposition 
2.3.9) 



\l i~f*i~f* f ry> ry* ^\ f ry ^ I ry*^ 

E u (t,x,x) = —z—Iv — exp 



(5.2) 



2t \2t J r \ At 

x u+i/2 eX p(-x 2 /4i) 



NE v (t,x) : = lim (x~ u ~ 1/2 E u (x, x, t)) 



where we have used the asymptotic behaviour of the modified Bessel function of first 
kind, cf. [AbSt92] 

(r/2Y 

'» {r) ~ WTT)' as r °- 

We define 

•t 



(5.3) 



F?(h)(t,x,y):=c v - [ [ NE v (t,x)H 1 g>(t,y,y)h(t - t,y) dtdy, 

JO JR b 



(-1), for u = 0, 
2u, for v G (0,1). 



for any h G L 1 (M + x M. b ) which is L^JO, oo) fl C°°(0, oo) in the time variable. Here 
H K b denotes the Euclidean heat kernel on M. b . The asymptotic behavior of F^(h) as 
x — > is studied by means of the Laplace transform Jzf in the time variable t G R + and 
the Fourier transform & in the Euclidean variable y G M 6 . For any g G L 1 (M + x M 6 ) 
which is ivi oc [0, oo) fl C°°(0, oo) in the time variable t G 1R + , not growing exponentially 
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as t — > oo, both transforms are defined as follows 

(5.4) {S?g){(,y)= f g(t, y) exp(-(t)dt, Re(C) > 0, 

</R+ 

(5.5) (&g)(t,u)= [ g{t,y)e-^dy. 

Jm. b 

Hence we assume henceforth that h e L 1 (M+ x R b ) and L^JO, oo) n C°°(0, oo) in the 
time variable t e R + , not growing exponentially as t — > oo. Finally, the inverse Laplace 
transform is given for any any 5 > and analytic L((), integrable over Re(£) = 5, by 

(5.6) (ST 1 L)(t) = ± f e*L(C)dC- 

Proposition 5.1. F^(h) is indeed the signaling solution to (5.1). In particular 

c-(F u N (h)) = h, c + (F„ N (h)) = G?*h(t):= f I G»(t-t,y-y)h(t,y)dtdy, 

JO J MP 

where (_? o &G N V ) (C, w) = C7^(C + M 2 ) 



(5.7) 



log VC+ |w| 2 + 7-log2, i/ = 0, 

S^2-^(c+M 2 r, ^(o,i). 

r i/ 



Proof. We compute for Re(£) > 

(JSf o &F?{h))(x, u, () = c u - %{NE V ■ &H Rb )(x, u, () ■ {J? o &h)((, u>) 

NE„(t, x)e~ t{ -^ 2) dt ■ o &h){(, u) 



o 



where i^T is the modified Bessel function of second kind, and in the definition of \/£ we 
fix the branch of logarithm in C\M _ . We have assumed that h(t, y) is not of exponential 
growth as t — > oo, so that (j£f o JPh)((,cu) is well-defined for Re(C) > 0. The Bessel 
function K v (z) admits an asymptotic expansion, see [AbSt92] 



K v {z) = {^,,^ _ v , _,_J _ n }+K(z),K(z)=O(z)asz^0. 



(log 2 — 7) — log(;z), for v — 0, 
2 v - 1 T{v)z- u + 2- y - 1 T{-v)z v , u>0, 

where 7 6 R is the Euler constant. Consequently 

(log 2 -log v /cTMF-7), ^ = 0, 
c + (^ o C, w)) = • (JSf o J^)(C, w) • <( r(-i/)(C + M 2 ) 



i> + 1) 2 2v+l 



(0,1), 



c-(j2f o J^f (/*)(•, C, w)) = o ^)(C, w). 
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Taking the inverse Laplace and Fourier transform, we obtain 
F u N (h)(x, t) = <fi{z)(Sf o &)-\ c + {Sf o &F?{h))) 

+ <j>z{x){se o &y\c-{se o &F?(h))) 

+ 2^+1) ^ ° ° * h ) K(xVcTW)(c + M 2 r /2 ), 

where each cf)^ coefficient exists and the third summand is 0(x 3 ' 2 ), 0. Conse- 

quently, indeed 

JSf o &(c±{F?{h))) =c±(ifc &F™{h)). 
This yields the stated explicit expression for the Laplace-Fourier transform of 



(5.8) (J?oJ?G»)((,u)={ V (- Ul , 



log VC+ |w| 2 + 7-log2, i/ = 0, 
(C + M 2 )", z^G (0,1). 



r(i/) 

Finally, we should note, evaluating the Fourier and the Laplace transform integrals 
explicitly, the following relation 

(5.9) G»(t, y-y) = H Rb (t iy - y)^ x G N u {t). 

□ 

6. Solution to the signaling problem 

We expand the lifted heat kernel (3*H& (as before rescaled under $ p from §2) asymp- 
totically at the left boundary face, using projective coordinates (4.3), where s = x/x is 
the defining function of If and x the defining function of the front face. We obtain by 
Theorem 4.2 

v>0 feeN 

where the sum runs over v > with v 2 being the eigenvalues of the tangential operator 
A in (2.3), and natural numbers k > 0. The coefficient G k v is a polyhomogeneous 
function on the left boundary face, of leading order (—1 — 6) at the front face and 
vanishing to infinite order at tf. In the projective coordinates (4.3), we find 

P*(xd x ) = xd x - sd s , f3*(xd y ) = d u , p*(d z ) = d z , f3*(x 2 d t ) = d T . 

Hence the action of f3*(x 2 (dt + A p )) keeps the order of s invariant. Consequently, since 
if jr solves the heat equation, so does each coefficient G k v s u+l l 2+k in the heat kernel 
expansion at If. We define 

(6.1) G° u s v+1 ' 2 = G° u F +1 / 2 x—V2 =: Hv F +i/2 ; 

where H v solves heat equation, since so does G° s u+1 ^ 2 . The kernel H u lifts to a polyho- 
mogeneous function on If of leading order (—3/2 — b — v) at the front face and vanishing 
to infinite order at tf. The left face If is a parabolic blowup of x M^ yjZ ) x 9M(^ at 

the highest codimension corner {t = 0, x = 0, y = y} and we denote the corresponding 
blowdown map again by 0. 
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Globally, H v is thus denned by the regularized limit in the sense that divergent 
summands are neglected 

P*H V = reg-lim(/3*F^ • Pl 7~ 1/2 ) ■ p ff ^ 1/2 . 

The relation between the kernels (5.2) in the model situation and the kernels H t ^, H v 
in the setup of an admissible edge manifold is then as follows. According to the heat 
kernel construction in §4, the expansion of the lifts P*H^ and (3*H U at ff in projective 
coordinates (4.2) at the top corner of * s given by 

/3*H*(p, £, I u, y,z,z) = p- 1 -' E »0-, f , l)P»(*> ^)H Rb (l, u) + 0*Kt + f3*K 2 
(6.2) »>o 

f3*H v (p } £, u, y } z } 7)= (p- 3 / 2 - b -»NE u (l, (1, u) + p* Kl + (5*k 2 ) P u (z } 5), 

as p — > 0, where if li2 and k 12 are the higher order terms, the sum runs over v > 
with v 2 being the eigenvalues of the tangential operator A in (2.3), and P u (z,z) is 
the Schwartz kernel of the fibrewise projection onto the corresponding eigenspaces. 
Moreover, as (£, £, p) — > 0, the higher order terms satisfy 



v>0 j=0 

(6.3) 



ex. 



(3*k 2 = J2J2°(p~ b+1 tt® u+1/2+i ^ p*** = o( P i/2 - b - u c +1/2 )- 

u>0 j=0 

We define for any v E [0,1) and h e l4 c ([0,oo) x 9M) n C7°°((0,oo) x dM), in 
standard coordinates in the collar neighborhood % C M 

F u (h)(t, x, y, z) := c v [ [ H u (t, x, y, z, y, z) h(t -t,y,z) dtdvol gM {y, z) 

=: c v H u * h. 



JdM 



Since H u solves the heat equation, so does F u (h). 

The fundamental component in the heat kernel construction of Mooers in [Moo99] 
is a solution u(t, •) 6 ^ max (A p ) to the signaling problem 

{d t + A p )u = 0, w(0,-) = 0, 

(6.4) 

where P u is the fibrewise projection onto the z/ 2 -eigenspace of the tangential operator 
A in (2.3). Note that for P u h(t) ^ 0, the solution u(t, ■) cannot lie in any self-adjoint 
domain of A p since by uniqueness of solutions to the heat equation u(0, •) = then 
implies u = 0. The signaling solution is in fact also unique, since for c~(u(t, •)) = 0, 
u(t, •) G @(Ap) and hence u(0, •) = then implies u = 0. 
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Theorem 6.1. F u (h),u G [0,1), is a signaling solution to (6.4). More precisely 

F v (h)(x) ~ c + (F v (h))tf{x) + c-{F u {h))^-{x) + 0(x 3/2 ), x -> 0, 

(6.5) ^(^(/i)) = [(C + G'J * P„/i + G" * + ^ • P„h] (t, y, z), 
c-{F v {h)) = P v h{t,y,z), 

where was introduced in Proposition 5.1 and G' U ,G" lift to polyhomogeneous func- 
tions on the parabolic blowup ofM> + x B 2 around Y := {(t, y,y) G M + x B 2 \ t = 0, y = y} 
of leading order (—2u—b) at the front face. 9 V G R is a constant and zero unless v = 1/2. 
The projective coordinates on [M + x B 2 , Y] near its front face are 

(6.6) P=Vi, u = — ;=-, y. 

Proof. The fundamental argument below is to explain how the higher order terms in the 
front face expansion of /3*H V contribute to the c ± (F u (h)) coefficients in the asymptotic 
expansion of F u (h) as x — > 0. Integrating first in z along the fibres, we obtain 

F v (h)(t, x, y,z) = c u I I NE u (t, x)H Rb (t, y - y)P u h(t - t, y, z) dtdvo\ B (y) 
Jo Jb 

+ c u / Kx{t,x,y,y^P v h(t-t,y,z)dtdvo\ B ($) 
Jo Jb 

+ c u / K 2 (t,x,y,y)P u h(t -t,y,z) dtdvo\ B (y). 
Jo Jb 

Expanding dvo\ B (y) in the first integral around y we find 

F v (h)(t, x, y,z) = c„ / / NE u (t, x)H Rb (t, y - y)P v h(t - t, y, z) dtdy 
Jo Jb 

+ c u I I (k + «i + K 2 )(t,x,y,y)P v h(t -%y, z) dtdvo\ B (y) 
Jo Jb 

=:F?{h) + F + F 1 + F 2 , 

where /3*k ,i = 0{ P ~ v ~ b ~ l/2 p v J~ 1/2 pg) and (3*k 2 = 0(p s u ~ b+1/2 p P + 1/2 p$) when lifted to 
If as a parabolic blowup of M + x M x DM at the highest codimension corner. 

The contribution to c ± (F u (h)) coming from the first integral F^(h) is studied in 
Proposition 5.1. It remains to discuss the latter three integrals -Po,i,2- Since we expand 
as x — > 0, for fixed t > 0, we may assume x 2 < t and separate for each j = 0,1,2 



,:2 



F 3 = j • J : P] • 
We describe the integral expressions Fj in the projective coordinates 

t y-y 

(6.7) r = — , u = , x, y, z, z, 

X z X 
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near the the lower corner of the left face If, where x = pg is the front face defining 
function and r = p t f the defining function of the temporal diagonal. Then, writing 
dvol b(v) = v(y)dy in local coordinates, we find 

Fj — / G'j(x,T,u,y, z)P v h{t — x 2 r,y — xu, z)v(y — xu) dr du 

JO JR b 



X 



x-" +3/2 , J = 0,1, 
x-^l\ j = 2, 



where G'j is bounded, polyhomogeneous and vanishing to infinite order as r — > 0, \u\ — > 
oo. Expanding P u h and v in Taylorseries around (t,y,z), as well as expanding G'j 
as x — > 0, we find no contribution to terms x ±v+l l 2 and y/x\og(x) in the asymptotic 
expansion as x — > 0, unless v = 1/2. 

In case v = 1/2, we have (— v + 3/2) = v +1/2 and hence Fj contributes • P^/i to 
the coefficient c + (F u (h)) if j = 0, 1, with no contribution to c~(F u (h)). In case v ^ 1/2, 
neither of Fj contributes to the coefficients ^(F^h)) and we set 6^- = 0. 

For the analysis of F",j = 0, 1, 2, consider the projective coordinates 

f- ^ x y — y ~ 

p=Vt, £ = -, u = , y, z, z, 

P P 

near the upper corner of the left face If, where p = pg defines the front face in these 
coordinates and £ = p r f the right boundary face. 

Contribution from Fq . By construction we find in the new projective coordinates 

Vt 

OO 

F o =E^[ / / NE u (l,OH Rb (l,u)p~^ 2 (pu) k v( k \y)PMt-p 2 ,y-pu,z)dpdu. 

k=l x MP 

Note that uH R b(l,u) = ^d u H R b(l,u) and hence integrating by parts in u, we find 

oo ^ 

F o =E^r// ^(i,0^(i,«)p-^ 1/2+ *t; ( * ) (i/)« fc - 1 



fe=i 



x (d y P u h)(t — p 2 ,y — pu, z) dpdu 

In view of the explicit structure NE V {1,£) = C £ u+1/2 = C ' x v+1 ' 2 p' v ~ x l 2 for some 
explicit constant C G R, we obtain 

F{ = x^ 2 J J p- 2 » +1 GZ(p,u,y,z)(d y P u h)(t - p 2 ,y - pu,z) dpdu = J J - J |, 

using the fact that p~ 2jy+1 is integrable at zero for v G [0, 1). Here, Gq is bounded in 
its components and polyhomogeneous in p. Obviously, the latter summand does not 
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contribute to the coefficients of x ±u+1 ^ 2 and y/x log(x) and hence we finally obtain (we 
abuse the notation by incorporating the p— factors into the kernel Gq) 

c-(F£) = 0, c + (F' { ;) = Gl* t d y P v h= [ [ G'^t,y,f)d y PMt-t,y,z)dtdvol B (y), 

JO J B 

where G' ' lifts to a polyhomogeneous function on the parabolic blowup of R + x B 2 
around Y := {(t, y,y) e IR + x B 2 \ t = 0, y — y} of leading order (— 2v — b) at the front 
face. 

Contribution from F". Here a more detailed information on the structure of k,i 
is necessary. Recall (4.8), which asserts that the second order term in the front face 
expansion of is given by J + * where J G ^f(M) with S> = 

(Ex, E ri + 1), and e *J5f (M) with 8 = (E lh E rf ). Consequently we may write 

Kl = kj + H {1) * d y K L , 

where kj lifts to a polyhomogeneous function of If of leading order (—1/2 — v — b) at 
the front face and of order (u + 3/2) at the right boundary face. Similarly, kl lifts to 
a polyhomogeneous function of If of leading order (—3/2 — v — b) at the front face and 
of order (y + 1/2) at the right boundary face. Consequently we may write 

F" = c u [ [ Kj(p,£,u,y)P v h(t - p 2 ,y - pu, z) p 1+b dpdu 

Jx JR b 

+ c u [ I H {1) * (p~ 1 d u + dy)K L (t - p 2 ,x,p,u,y) P v h(t - p 2 ,y - pu,z)p l+h dpdu, 

Jx JR b 

where we have neglected the factor v(y) in the volume form dvol#(y) = v(y)dy. Inte- 
grating by parts in the last integral we arrive by the composition law in Theorem 4.4 
at the following expression 

F'/ = x^' 2 [ 1 [ p- 2 ^Gj(p,i,u,y)P v h{t- p 2 ,y- pu,z)dpdu 

Jx JR b 

+ x v+1/2 I [ p- 2v+l G L , a (p, £, u, y)P v hit - P 2 ,y~ pu, z) dp du 
Jx Jr>> 

+ x» +1 ' 2 [ [ p- 2u+l G L ^p,i,u,y){d y P v h)(t - p 2 ,y - pu,z)dpdu, 
Jx Jr>> 

where the kernels Gj,Gl^,Gl^ are bounded and polyhomogeneous in p. 

For the first integral, expanding P u h in Taylorseries around (t,y,z), as well as ex- 
panding Gj as p — > 0, we find no contribution to terms x ±u+l ^ 2 and y/x\og{x) in 
the asymptotic expansion as x — > 0, unless v — 1/2. In case v = 1/2, we have 
{—v + 3/2) = v + 1/2 and hence the first integral contributes 9" ■ P v h to the coefficient 
c + (F v (h)), with no contribution to c~(F u (h)). In case v ^ 1/2, the first integral does 
not contribute to the coefficients c^^F^h)) and we set 9" = 0. 
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For the latter two integrals the discussion is parallel to that of Fq and we obtain (we 
again abuse the notation by incorporating the p— factors into the kernels G>) 

c~(F") = 0, c + (F") = G Lia * P u h + G L , 6 * d y P u h + const ■ P v h 

= G L>a (t,y,y) P v h(t -t,y,z)dtdvo\ B (y) 

Jo Jb 

+ / / G Ljb (t,y,y)d y P u h(t-t,y,z)dtdvol B (y) 

Jo Jb 

+ 6'{-P„h{t,y,z) 

where Gi >a , Glj, lift to a polyhomogeneous function on the parabolic blowup of 1R + x B 2 
around Y := {(£, y,y) G IR + x B 2 \ t = 0, y = y] of leading order (— 2u — b) at the front 
face. 

Contribution from F%. As before we obtain by construction 

F' 2 = c v \ \ K 2 (t,x,y,y) P v h{t -t,y, z) dtdwo\ B (y) 
Jx 2 Jb 

= x» +1 ' 2 [ * / p- 2v+1 G'l(p,Z,u,y)P v h(t- P\v~ pu,z)v(y- pu)dpdu 

Jx JR b 



' x 
111 • 



where the kernel G 2 is bounded and polyhomogeneous in p. Its discussion is parallel to 
that of Fq and we obtain (we again abuse the notation by incorporating the p— factors 
into the kernel G" 2 ') 

c-(F^) = 0, c + (F») = G'i * P v h = f [ G»(t, y, y) P v h(t - t, y, z) dt dvo\ B (y). 

Jo Jb 

where G 2 ' lifts to a polyhomogeneous function on the parabolic blowup of M + x B 2 
around Y := {(t, y, y) G R + x B 2 \ t — 0, y — y} of leading order (— 2u — b) at the 
front face. This proves the statement with G' v = Gi a + G^, G" = Gq + G^b and 
0*, = O + 0' x + 0?. D 



7. Heat kernel for algebraic boundary conditions 



In this section we finally employ the signaling solution to construct the heat kernel 
for a given algebraic boundary conditions T. Consider the increasing sequence of eigen- 
values v 2 G [0, 1) with j = 1, ...,p, of the tangential operator A in (2.3), counted with 
their multiplicities. Consider G C™(M) and put u = Hj?(f). We seek to correct u to 
satisfy algebraic boundary conditions T, i.e. 

v 

(7.1) w = u + ^F Vj {h Vj ) G^(Ar). 

j'=i 

where each h u . lies in ImP^,.. Recall T = (Ty) G Matr(g, A q ) with diagonal entries 
given by Tjj = bjjipj +0jjipj~, and the off-diagonal entries = The coefficients 

bij,0ij G K. are such that either ba = 1, or ba = 0, where in the latter case we require 
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On = 1 and 8ij = for i ^ j. Here we assume that bjj = 1 for every j = 1, ..,p, since in 
case of bjj = 0, h Vj = c£ (u). Then (7.1) reads as follows 

p 

(7.2) c+(«) +c + (F^(^.)) = J^A,, i = 

i=i 

We define p x p matrix valued operators 

G w := -diag(G^,...,G^) - diag(^, ... AJ. 

g : = diag(G / , i + g^, g; + g;<9,). 

We also write 

H v 4> := diag(if yi 0, .., #„ p 0), 
h v := diag^, ..,/i^). 

Then we may rewrite (7.2) as H u <p = (G N — G)h u . We apply the Fourier transform in 
y G M. b and the Laplace transform in t G M + , and obtain 

if o &H v <j> = %o &G N ■ if o &h v -if o ^(G?^) 

(7.3) = if o ^G w • (if o - (if o ^G^)- 1 • if o &{Gh v )) 

= G N ■ o - G^ 1 • if o ^(G/i, 

where we have set 

G N := o = (%)„ - diag(i^ o J?G^, if o ^G^) - diag(^ 1 , 9 V 



Recall from Proposition 5.1 that Gn is a function of (C + M 2 )- Existence and the 
structure of the inverse matrix Gn(( + M 2 )" 1 is the core of the subsequent section. 



Put 



D := o if- 1 G N X if o 



Then, applying the inverse Gn(( + \w\ 2 ) 1 on both sides of (7.3), we arrive at the 
following relation 

(7.4) DH v (j) = h u — D o Gh u . 

In the setup of isolated conical singularities, the operator G is absent and (7.4) provides 
an explicit result for h u , obtained in one step by inverting the corresponding matrix. 
Unlike in the case of isolated conical singularities, on edges G is a non trivial operator on 
the base manifold B and the solution is obtained from (7.4) by an iterative procedure. 
We define for any M G N 

^ m \ r m -j 

^(DoGfoD H u 

_k=0 



'Vi 



(7.5) 
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This defines an approximate solution to (7.4), solving it up to an error 

(7.6) DH v <f> = h A J -Do Gh A J -(Do G) M+1 DH u <f). 

Below we show that the error term vanishes and h\f converges as M — > oo. Anticipating 
that discussion, the heat kernel for algebraic boundary conditions T is subsequently 
given by 

p 

(7.7) H r = H^+J2 H »i * ( K rh * H »i , 

where the kernels are convolved in t and concatenated over B. 

The remainder of the section is concerned with the analysis of the operator orders in 
the definition of h M . The integral kernel of D is explicitly given by (u = (y — y)/yt) 

D(t,y,y)= [ [ e*b-®tXG£(C + \ U \ 2 )dCdu 

(7.8) =\\ e- iu{y -y ) e- t ^ 2 e st G- N 1 {s)dsdu 



The asymptotics of Jz? 1 G N 1 (t) as t — > is established in Theorem 8.2 below. More 
precisely, noting the p x p matrix structure of the operators, we obtain for any (ij) G 
{1,-,P} 2 

where Uij — Vi + Vj if i ^ j and Vij — v t if % — j. Consequently, cf. the projective coor- 
dinates (6.6), the component D^ lifts to a polyhomogeneous function on the parabolic 
blowup of R + x B 2 around Y := {(t,y,y) G M + x B 2 \ t = 0, y = y} of leading order 
(—2 — b + 2uij) at the front face. 

Proposition 7.1. Let Gi and G 2 lift to polyhomogeneous functions on the parabolic 
blowup o/R + x B 2 around Y := {{t,y,y) G 1R + x B 2 \ t = 0, y = y} of leading orders 
(—2 — b + «i) and (—2 — b + a 2 ) at the front face, respectively. Then 

G 1 * t G 2 (t,y,y') = / G^t -t,y,y)G 2 (t,y,y') dtdvo\ B {y), 

JO J B 

lifts to a polyhomogeneous function on [M.~^ x B 2 ,Y] of leading order (—2 — b + ai + a 2 ) . 

Proof. By polyhomogeneity of we may write for i — 1,2 and any N G N, c > 

( G*(c 2 t,cy,cy) = C - 2 - b ^ +k G^t,y,y), 
G?(c 2 t,cy,cy) = 0(c- 2 - b ^+ N ), c -+ 0. 



Gi(t, y,y) = J2 ° k i& v>v) + G f (*> y>y)> 



k=0 
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For a composition of individual homogeneous summands we obtain 

G\ * t G l 2 (c\ cy, cy') =f'f G\ (cH - t, cy, y) G l 2 (t, y, cy') dtdvo\ B (y) 

JO J B 

= c 2+b [ [ G k (c 2 t~c 2 t,cy,cy)G l 2 (c 2 t,cy,cy')dtdvol B (y) 
Jo Jb 

= c- 2 - b+a ^ +k+l G k 1 *tG l 2 (t,y,y'). 

Similar estimates for the remainder terms G± 2 prove the statement. □ 

The statement on the leading orders in Proposition 7.1 above holds also without the 
assumption of polyhomogeneity. Moreover, presence of d y derivatives is not excluded 
from that picture, since d y corresponds to lowering the leading order at the front face of 
[M. + x B 2 , Y] by one. Consequently, G is zero off diagonal and its diagonal components 
Gjj = (G' v . + G" v .dy) may be viewed as of leading order (—1 — b + 2Vj) at the front face, 
when lifted to [R+ x B 2 ,Y]. 

Hence, by Proposition 7.1, we find that {D o G)ij = Dij o Gjj lifts to a poly- 
homogeneous distribution on [M + x B 2 ,Y] of leading order (—2 — b + where 
otij := 1 + 2(z/y — Vj) > 1. Due to > 1, subsequent M-fold compositions improve 
the order and hence h^f converges and the error in (7.6) vanishes as M — > oo. 



8. Heat-trace asymptotic expansions 

In view of (7.8), we begin with studying G^ 1 and its inverse Laplace transform. Fix 
the branch of logarithm in C\1R~ for the definition of powers of complex numbers. We 
apply the inversion rule using adjuncts 

(8.i) (d-^j = (-iy+ > djG ^ (c) , 

y } detG JV (C) 

where (adjGiv)ij(C) is the determinant of the reduced matrix, obtained from Gn by 
deleting the z-th row and j-th column. Set Kg := 2(7 — log 2 + 6) and introduce the 
multi- index notation, setting for a = (a±, a q ) G and ft = (ft q +i, ■■, ftp) G W~ q 

(iogc+^r = noogc+K*j afc , ^= n c^ fc . 

k=l k=q+l 

Let tuples with each entry given by 1, be denoted by 1. Then in view of the explicit 
formulas in Proposition 5.1 we obtain 

(q p 
l«l=o |^|=o 
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for certain coefficients C,C a p. The summation above excludes \a\ = = 0. For 
\(\ 3> we may expand (det Gjv(C)) -1 m Neumann series and obtain 



(detGW(C)) 



-i 



cr^OogC + ^Cv)" 1 [i+ E A*(iogC + **)~ a C /s 



which is convergent for |£| 3> and in particular asymptotic series as \(\ — > oo. Similar 
computations hold for (&d]F)ij((). Hence overall we obtain, cf. ([Moo96], Lemma 8.6) 

Proposition 8.1. There exist constants Aij and B a p such that for \(\ ^> 



V M,l/3|=i 

(logC + ^) _1 ( lo gC + if v% = = 0,i ^j, 

(logC + ^)- 1 C" !/j , if v i = 0,v j jL0, 



x < 



(logC + ^i)- 1 , 



if z/i 



0, i = j, 



if v u Uj ^0,i^ j, 
if Vi = v i ^ 0, % = j. 



We point out that no (log £) factors arise if we have Friedrichs extension on the v = 
components of the Laplacian. This case exzibits only a classical expansion without the 
exotic phenomena and we do not consider it here. We now want to derive an asymptotic 
expansion for (5 > fixed) 

2>- x {x,m = ^-.( e K Gogc + ^r a c^c 

Z7n JiR+S 

Theorem 8.2. Write v$ := P q+ \V q+ i + .. + f$ v v v for the given (3 G W~ q . For some 
a G N 9 and j3 G N p_g we then have up to smooth additive components 

oo 

~ X>0 fc r 1+ "'io g -w-*(*), t -> o. 

fc=0 

Proof. As the first step we deform the integration region (zR + 5) to //, concatenated 
out of three parts, /ii = z(— oo, 1], the half circle /i 2 = {C G C | |£| = l,Re(£) > 0} 
oriented counterclockwise, and /X3 = z[l,oo). The change of the integration contour is 
possible, since 

r-5 



2m J 



< (}ogR)- lal R-^ -»■ 0, as i? 00. 



Integration over /j 2 leads to a function 1 (^ i ^) At2 G C°°[0, 00), so that we may write 

1 



2vr« 
1 



2tt 



MlUM3 



2tt 
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For the second summand we deform the integration contour to oo). For the third 
summand we deform the integration contour to (— oo)). We explicate the argument 
for the second summand. Let R > 1 and the contour r/ R := {i?exp(i0) | G [0, 7r/2]} 
be oriented counterclockwise. 





^\ Vr 




R 



Figure 2. The integration contour rj R . 

We use the O-notation for the asymptotics as R — > oo. Then, substituting x 
Rexp(i(p), we find for some C > 

/ I e itx XA ix ) | dx = R I 1 | exp (-tR sin (j) + i((f) + tR cos ^^(iRe^) \ dcf. 

Jn„ Jo 



<CR(\ogR)-^R-^ / exp(-tR sin <j>) dxf> 

Jo 

= C R (log R)~ H R^ / exp(-ti2sin0)# + O(ir°°) 

Jo 



/•tt/4 

< C R (log R)~^R~^ I exp(-tR sin 0) cos(0) d(j) + 0(iT°°) 

Jo 

psimr/4 

= C R (log Ry^R-^ / exp( y -tRy)dy + 0(R- co ) 

Jo 

= O ((log Ry H R~^) -> 0, as R ->■ oo, 

where in the fourth line we have used the fact that cos</> is bounded from below for 
<f) G [0,7r/4], and in the fifth line we substituted y = sin0. Subsequently, writing 
fji = {z | z G for the clockwise oriented contour, we find 



-Uar \u\-cs-Uar ■ ' I e itx XA ix ) dx 



H2 



2tt 



m 



e- ltx £r aP (-i x )dx 



+ e- t y^ a A-y) + dy--j^ e^^A^-dy, 

where (— y)± denotes negative real numbers with the argument (±7r), and the first three 
summands define smooth functions C°°[0, oo). 
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We continue under equivalence up to smooth functions. Then without explicating 
the smooth summands we obtain 



ZTT 



I 

2n~ 



J2 C i [°° e-^y-^ilogy + tn + Ke^dy) 

/oo 



We establish an asymptotic expansion for each Jzf^y(t) as t — > 0. For this we differ- 
entiate [up] times in t and obtain by a change of variables x = ty 

-^K^) = C % / e" V^Oogy ± in + Ke)- p dy 

poo 

= ^(-1)^ ] Cf p / e- ty y- v e +[v t> ] (\ogy ± in + K e )~ p dy + smooth 

x r 1+ ^~ [ ^ ] log" |p| (t) + smooth 

Since (1 — r) _1 = X]fclo rfc + r M+l {l — r) _1 for any M E N and r ^ 1, we find, up to 
smooth additive components 

~ r 1+ »'-w iog" |a| (t) ^^ fc iog- fc (t), t -> o. 

fc=0 

In order to integrate the asymptotic series [i/g] times, note for any M > — 1 and iVeN 

fr M \og- N {r) dr = f{M + 1)~ V' +1 \og~ N {t))' dr + jT ^^^"V) <*r 

= (M + 1)-4 M+1 bg-^(t) + 0(t M+1 log-^ 1 ^)), t 0. 

Iterating this argument we finally arrive, up to smooth additive components, at the 
asymptotic expansion 

oo 

s#(t) ~r 1+ "'b g -W(t) io g - fe (t), *->o. 

fc=0 

This proves the statement. □ 

Let K lift to polyhomogeneous function on the parabolic blowup of M + x B 2 around 
Y := {(£, y,y) e R + xB 2 | t = 0, y — y} of leading order (—2—6+7). By an ad verbatim 
extension of Proposition 7.1, compare also ([Moo96], Proposition 8.7), we find that 
the convolution H Vi * K * H v lifts to a polyhomogeneous distribution on the parabolic 
blowup of M+ x (dM) 2 around Y F := {(t, (y,z) } (y,T>) G M+ x (dM) 2 \t = 0,y = y}, 
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of leading order (— 1 — b — (z/j + Vj) + 7) at the front face, and vanishes to infinite order 
at tf and td. By standard pushforward arguments we obtain 



5.2) 



Tr H Vi *K* H v . (t) ~ Vi 7 ^ d k Vi \ as t ->■ 0. 



fc=0 



Combining Theorem 8.2 with (8.2), we obtain the following 

Proposition 8.3. Let K lift to polyhomogeneous function on the parabolic blowup of 
W + x B 2 around Y := {(t,y,y) G 1R + x B 2 \ t = 0,y = y} of leading order (—2 — 6 + 7). 
Write up := (3 q+ iv q+1 + .. + (3 p v p for the given (3 G W~ q . For any a eW and (5 G W~ q 
we put So.fi = (logC + Ke)~ a (v ■ Then as t — > we obtain 



Tr [H Ui * ST\^ P ) * K * (t) ~ Vi 



J2d k V~t k \og-^-\t). 



k,l=0 



Proof. The statement follows by an iterative application of the following identity 

j\«\Og'(T)(t - tYAt = £(-1)* Q j\°+ k \0Z>(T)dT 



E 

k=0 

00 

E 

fc=0 



(~l) fc (ji. 
a + k + l\k 

(-D h 



1 (fc) l0g _P(t) + ^S"^ 1 ^))) ' t °- 



□ 



Combining Proposition 8.1 with Proposition 8.3 we obtain a full asymptotic expan- 
sion for each TrH Ui * t (K^)ij *tH V] , complicated only by the intricate combination of the 
various components after matrix multiplications. The heat trace asymptotic expansion 
admits logarithmic terms in accordance to [KLP08] and [GKM10]. We explify the 
leading orders in our final main result. 



Theorem 8.4. 



f log- 2 (t), 



if Vi 



Vi 'log" 1 ®, if 1/4 = 0,1/^0, 



TrH Vi * t * t H v . ^ { 



log-\t), 

Vi , 



if Vi 



V; 



0, i = j, 



if v^ uj ^0,i^ j, 
if Vi = uj ^ 0, i = j. 



This corresponds to ([Moo96], Theorem 8.2), up to the case z/j — Vj 7^ 0, i = j. 
Note for example that for z/j, Uj 7^ and i ^ j the correcting kernel lies according 
to ([Moo96], Theorem 8.2) in the space $"«+ I 'i+ 2 ) with the front face lading order 
—3 + Vi + Vj + 2 = z/j + z/j — 1, which leads to f^+^O/ 2 after taking traces. Here, 
Theorem 8.4 and ([Moo96], Theorem 8.2) agree. 
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However, in case i = j, asymptotics of G^ 1 ^) m Proposition 8.1 is different and hence 
Theorem 8.4 states t° expansion, different from the expansion for i ^ j, in contrast to 
([Moo96], Theorem 8.2) which still asserts a heat trace expansion of order t^+^j)/ 2 . 
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